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Abstract 

We apply the gradient approach to obtain a path over the sphaleron barrier and 
to demonstrate the fermionic level crossing phenomenon. Neglecting the mixing 
angle dependence and assuming that the fermions of a doublet are degenerate in 
mass we employ spherically symmetric ansatze for the fields. The gradient path 
over the barrier is smooth, even for large values of the Higgs boson mass or of the 
fermion mass, where the extremal energy path bifurcates. 
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1 Introduction 



In 1976 't Hooft observed that the standard model does not absolutely conserve baryon 
and lepton number due to the Adler-Bell-Jackiw anomaly. The process 't Hooft considered 
was spontaneous fermion number violation due to instanton induced transitions. Fermion 
number violating tunneling transitions between topologically distinct vacua might indeed 
be observable at high energies at future accelerators 0, |3|. 

Manton considered the possibility of fermion number violation in the standard model 
from another point of view ||. Investigating the topological structure of the configuration 
space of the Weinberg-Salam theory, Manton showed that there are noncontractible loops 
in configuration space, and predicted the existence of a static, unstable solution of the field 
equations, a sphaleron || , representing the top of the energy barrier between topologically 
distinct vacua. 

At finite temperature this energy barrier between topologically distinct vacua can be 
overcome due to thermal fluctuations of the fields, and fermion number violating vacuum 
to vacuum transitions involving changes of baryon and lepton number can occur. The rate 
for such baryon number violating processes is largely determined by a Boltzmann factor, 
containing the height of the barrier at a given temperature and thus the energy of the 
sphaleron. Baryon number violation in the standard model due to such transitions over 
the barrier may be relevant for the generation of the baryon asymmetry of the universe 

While the barrier between topologically distinct vacua is traversed, the Chern-Simons 
number changes continuously from Nqs = in one vacuum sector to Ncs ± 1 in the neigh- 
bouring vacuum sectors, passing through the sphaleron at Ncs = ±| [0, 0- However, 
for large values of the Higgs boson mass, energetically lower, asymmetric sphaleron solu- 
tions appear, the bisphalerons ]T3| , 13 ]. The minimum energy path over the barrier | jTT| 
then develops bifurcations |l2j , indicating the need for another approach to the sphaleron 
barrier, which yields smooth paths. 

As the barrier is traversed one occupied fermion level crosses from the positive contin- 
uum to the negative continuum or vice versa, leading to the change in fermion number. 
When considered in the background field approximation this level crossing phenomenon 
predicts the existence of a fermion zero mode precisely at the top of the barrier, at the 
sphaleron [JIS, [16], [17]. For massless fermions this zero mode is known analytically |T5|, [Ll; 



Considering the minimum energy path over the barrier |TT], [L2[] the fermionic level 
crossing was demonstrated recently in the background field approximation under the as- 
sumption, that the fermions of a doublet are degenerate in mass [fl8 , 19 , |20[| . This as- 



sumption, violated in the standard model, allows for spherically symmetric ansatze for 
all of the fields, when the mixing angle dependence is neglected (which is an excellent 



approximation |21], p2|). An analogous, but selfconsistent calculation p3[ led to similar 
results. However, for heavy fermions it led to strongly deformed barriers, eventually giving 
rise to bifurcations and to new sphalerons . 



Motivated by the catastrophes encountered along the energy barrier in the extremal 
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energy path approach for large Higgs boson or fermion masses we here consider the gradient 
approach to the sphaleron barrier. In section 2 we briefly review the Weinberg-Salam 
Lagrangian and the anomalous currents for vanishing mixing angle and for degenerate 
fermion doublets. In section 3 we discuss the sphaleron barriers. We present the radially 
symmetric ansatz for the boson fields and obtain the energy functional. We discuss the 
gradient approach, putting special emphasis on the question of the underlying metric. We 
compare the barriers obtained with gradient approach with those of the extremal energy 
path approach. In section 4 we discuss the fermionic level crossing along the gradient path 
over the barrier. We present our conclusions in section 5. 



2 Weinberg-Salam Lagrangian 

We consider the bosonic sector of the Weinberg-Salam theory in the limit of vanishing 
mixing angle. In this limit the U(l) field decouples and can consistently be set to zero 

£b = -\F; v F^ a + (D^(D^) - A($t$ - I v 2 ) 2 (1) 
with the SU(2) £ field strength tensor 

F a , v = W - d v V; + ge ab XV u c , (2) 
and the covariant derivative for the Higgs field 

D^=(d»- l -i 9 T a vf)<S> . (3) 

The SU(2)l gauge symmetry is spontaneously broken due to the non-vanishing vacuum 
expectation value v of the Higgs field 



<•»-*(:)• 



(4) 



leading to the boson masses 



M W = M Z = l -gv , M H = vV2\ . (5) 

We employ the values M w = 80 GeV, g = 0.65. 

For vanishing mixing angle, considering only fermion doublets degenerate in mass, the 
fermion Lagrangian reads 

£f = qih^D^qL + q R i^d^q R 

- f {q) q L ^u R + M R ) - f {q \d R & + u R &)q L , (6) 
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where denotes the lefthanded doublet while qn abbreviates the righthanded 

singlets (ur,cIr), with covariant derivative 

D,q L =(d ll -Ugr a V^q L , (7) 

and with $ = ZT2$*. The fermion mass is given by 

M F = . (8) 

Due to the U(l) anomaly baryon number and lepton number are not conserved 

= -f^K, , (9) 

where 

K, = -|^ W Tr(^V CT + \igV u V p V°) (10) 

(F U p = l/2r l F* p , V ct = l/2r l V^) is the Chern-Simons current and / g is the number 
of generations. In the unitary gauge the topological baryon number Q B , carried by a 
configuration, is determined by its Chern-Simons number iVcs ; 

N cs = J d 3 rK° . (11) 

For the vacua the Chern-Simons number is identical to the integer winding number, while 
the sphaleron at the top of the barrier carries half integer Chern-Simons number ||. 



3 Barriers 

Approximations to the sphaleron barrier can be obtained by constructing families of field 
configurations for the gauge and Higgs boson fields, which interpolate smoothly from one 
vacuum sector to a neighbouring one. The minimum of all maximum energy configurations 
encountered along such vacuum to vacuum paths represents the sphaleron [Q or, for 
large Higgs boson masses, the bisphaleron |H| [TJJ. Applying the gradient approach, 
we construct the sphaleron barrier for various Higgs boson masses. We compare these 
barriers to the ones obtained by constructing the extremal energy path JIT] , 



3.1 Energy Functional 

In the limit of vanishing mixing angle the gauge and Higgs boson fields can be parametrized 
by a spherically symmetric ansatz given by 



v , = izM^ + MV-^ + M^., (12) 

gr gr gr 

Vq = > (13) 

$ = -^=(H(r)+ir-rK(r)) ( J J , (14) 
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which involves the five radial functions fA{ r ), f'B{ r ), fc{ r ), H(r) and K(r). 

This ansatz is form-invariant under spherically symmetric gauge transformations with 
the SU(2) matrix 



0(r) _ 
U(r) = exp(i f ■ r) 



(15) 



The functions then transform as 

fA + ifB — ► exp(iQ)(f A +if B ) , 
fc — fc + rQ', 

H + iK — ► exp(z ® ) (H + iK) . 
The ansatz, Eqs. fll2])-(Il), leads to the bosonic energy functional 



(16) 



9 A 



dx 



X 



) 2 + (f B -^f 



X 



fl 



where 



[K 2 + H 2 ){1 + fi + f* + i£) + 2f A {K 2 - H 2 ) - Af B HK 
2xf c (K'H - KH') + 2x 2 (H' 2 + K' 2 ) + ex 2 (H 2 + K 2 - if 



x = M w r 



(17) 



is a dimensionless coordinate, the prime denotes differentiation with respect to x, and 

(19) 



4A _ 1 / M H \2 
^2 - 2~V M^7/ ' 



w ' 



This energy functional has a residual U(l) gauge invariance with respect to the gauge 
transformation, Eq. (|i~6"l) . 

The Chern-Simons number of a given configuration is 



A, 



cs 



2vr Jo 



dx 



(fl+f 2 B)(--0') 



X 



■ fc 

' X 



&) 



{f\ + f 2 ) sin(# - 6) 



(20) 



where 

9(x) = arctan(/ B (x)/7 A (x)) . (21) 

The function Q(x) is an arbitrary radial function, associated with the U(l) gauge transfor- 
mation, Eq. (0). From the expression fl20|) the Chern-Simons number is readily obtained 
in an arbitrary gauge. 

In the radial gauge, where fc = 0, the spatial part of the Chern-Simons current 
contributes to the topological baryon number. One therefore has to rotate to the unitary 
gauge, where only the Chern-Simons number determines the topological baryon number. 
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The corresponding gauge transformation involves the function Q(x), which satisfies O(0) 
and ©(oo) = 8(00). This leads to the Chern- Simons number 



1 f°° , , 

Ncs = 7T dx (fBf' A ~ /a/b 
Ztt Jo 



00 



2tt 



:22) 



The energy functional, Eq. (17), possesses non-trivial extrema. The sphaleron [EJ 
with iVcs — 1/2 exists for all Higgs boson masses. For large Higgs boson masses the 



energetically lower, asymmetric bisphalerons bifurcate from the sphaleron |L3|, |14 . 



3.2 Gradient Approach 

Let us now consider the energy barriers, associated with the sphaleron and, for large 
Higgs boson masses, the bisphalerons. For instance, the minimum energy path over the 



sphaleron barrier is obtained from the functional [11 



W 



&b H n — 4^ v cs 



(23) 



where £ is a dimensionless Lagrange multiplier 
Higgs boson masses 
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This path is satisfactory for small 
For large Higgs boson masses, however, variation of the func- 
tional, Eq. (p3|), leads to an extremal energy path, which does not culminate at the 
bisphaleron, but has a spike in the vicinity of the bisphaleron and culminates at the 
sphaleron ||12|| . This is in clear contrast to expectation and to other paths constructed 
in more or less ad hoc fashion |2"0[ . Similarly, when the path is constructed in the 
presence of fermions with large masses, bifurcations arise along the path in the vicinity of 
the sphaleron |23[ . 

These catastrophic features of the extremal energy path are artifacts of this approach, 
indicating the need for another systematic approach to obtain the sphaleron barrier. 



3.2.1 Explicit gradient formalism 

Let us therefore consider the gradient approach as an alternative approach to the sphale- 
ron barrier. Starting at the sphaleron or bisphaleron we are looking for the steepest path, 
connecting the top of the barrier with the vacua on both sides. 

The direction of steepest descent at the top of the barrier is determined by the negative 
mode of the sphaleron or bisphaleron W\ |26|. Away from the top of the barrier, the 



gradient of the energy functional determines the direction of steepest descent. Thus for a 
given configuration / the neighbouring configuration / = f+Sf along the path of steepest 
descent is obtained by choosing 5f proportional to the gradient of the energy functional. 
Explicitly, when / denotes a set of functions fi, we obtain 5f from the functional derivative 
of the bosonic energy functional according to 



Sfi 



a 



d£(f,f) ( dS(f,f) 
dfi(x) { dfl{x) 



(24) 



/=/,/'=/' 
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with 



E h = J dx£(fi,fl) 



(25) 



and a is a small negative number. 
Metric 

The notion "steep" always refers to a metric on the configuration space. Therefore, 
to unambiguously define the gradient approach, we need to specify a metric on the con- 
figuration space. A natural choice for this metric is provided by the kinetic energy term 
of the Lagrangian [27 . For the spherically symmetric ansatz, Eqs. (^)-(JTJ), the effective 
mass 



28], [29], |30|, [31], [i2[ Q of the gauge-Higgs-system reads 

2 



m(A) 



•57T 



9 2 M, 



w 



dx 



dX 



'df B \ 2 , 1 / ~ 



dA 



2 \ dX 



dX 



dX 



(26) 

where A is an arbitrary path parameter. 

Accordingly we define a distance d of two configurations / = (f A , f B , fc, H, K) and 
/ — if A) /bj Jc> -H") -K"), taken at the "times" A and A, as 



rf 2 (/,7) 



16tt 



(/a - /a) 2 + Ub - /b) 2 + -{fc - fc) 



+ 2x 2 (F - Hf + 2a; 2 (ir - Rf 



(27) 



Since the gradient formalism assumes a Euclidian metric with equal weight for all indices 
(i. e. space points and function indices), corresponding to a distance d 



(28) 



x 



we have the relations f x {x) = f A (x), fo{x) = f B (x), / 3 (x) = ■j%fc(x), h(x) = V2xH( 

and /^(x) = \/2xK(x). 
Equations 

Assuming the old configuration / = (f A , f B , fc, H, K) is in the radial gauge, fc = 
the new configuration / = (f A , f B , f c , H, K) will in general not be in the radial gauge 
/c/0. Denoting 8f = f — f, we find the set of equations 



8f A = 2a ^(f A + f B -l)-f; + f A (K 2 + H 2 ) + K 2 -H> 
l x z 



f 



8f B = 2a ^(fl + f*-l)-f>>+f B (K* + H 2 



2HK 



c 



4a 



x 



f' A f B -f' B f A -x 2 {K'H-H'K) 



(29) 

(30) 
(31) 
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6H = 

a 



6K 



x' 2 



H(l +fl + f B ) - 2f A H - 2f B K + 2Hex 2 (H 2 + K - 1) - (2x 2 H')'\ , (32) 
K{1 + fl + fl) + 2f A K - 2f B H + 2Kex 2 (H 2 + K 2 - 1) - (2x 2 K')'] . (33) 



To obtain the new configuration in the radial gauge, / = (f A , /b, fc, H, K), we perform 
a gauge transformation, Eq. ([16]), of / = (f A ,f B ,fc,H,K) with the gauge function $ 
determined by /c = = fc + 

3.2.2 Gradient formalism with constraint 

Let us now consider a modified gradient formalism, which is equivalent to the explicit 
gradient approach in the limit d(f, /) — > 0. In the explicit gradient formalism the gradient 
is taken at the 'old' configuration /. In the modified gradient formalism we take the 
gradient at the 'new' configuration / = / + 5f. Then Eq. ( p4|) is a set of differential 
equations for the functions fi, which can also be obtained by variation of the functional 
W 

W(f) = E h (f) + ^M w d 2 (f, f) , (34) 

where a constraint is imposed on the distance d(f,f), and £ is a Lagrange multiplier 
inverse proportional to a. 
Choice of gauge 

As before, starting from a configuration in the radial gauge, fc = 0, the neighbouring 
configuration will in general not be in the radial gauge, /c / 0. To be able to keep the 
radial gauge throughout, we define a new distance d(f, f) 

d(f, f) = min d(f, /_$) = 



16tt 



mm 



i2 <j, 



J dx (f A - (/acos($) + /Bsin($))) 2 + (f B - (/bcos($) - /Asin($))) 2 



I 9 

+2x 2 (H - (i/cos($/2) + Ksin($/2))) 2 + 2x 2 (k - (iTcos(<l>/2) - sin($/2 
+^(fc-fc + x&) 2 



(35) 



where is obtained by gauge transforming / with the gauge function $, which minimizes 
the distance (^). In contrast with the distance, Eq. (|27|), the new distance is gauge in- 
variant under independent gauge transformations of / and /. (Denoting the configuration 
space by C, this is a proper metric on the projective space of gauge orbits CJ ~, where the 
equivalence relation ~ identifies configurations which are connected by radially symmetric 
gauge transformations, Eq. flT6"|) P3fl .) 
Equations 

Variation of the functional W, Eq. fl3"4]), with distance d(f, f) replaced by d(f, f) leads 
in the gauge fc = fc = ® f° r the gauge and Higgs field functions to the set of equations 

k 

X 1 



n= t 4Ul + fl-^) + fA{K 2 + H 2 )+K 2 -H 2 + i{f A -{f AC o^) + f B ^))) , (36) 
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f^ = ^(fl + f^-l) + fB(K 2 + H 2 )-2HK + ^(f B -(f B co S ^)-f A smm) , (37) 

Ob 

H" = - 2 -H' + ^- 2 ((l-f A r + f B )-^f B + ^ + K^-l)H 

+ t(H- (#cos($/2) + Xsin($/2))) , (38) 

K" = - 2 -K' + ^- 2 ((l + f A r + f B )-^f B + e(H^ + K^-l)K 

+ Z(K- (Xcos($/2) -ffsin(f/2))) , (39) 
and to an additional equation for $ 

= -^ , + ^sin($)(/ J 4/ A + / B / B ) + ^cos($)(/B/A-/A/i ? ) 

+2 sin($ /2) (M + M) + 2 cos($ /2) (M - HK) . (40) 

Boundary conditions 

The boundary conditions are chosen such that both the energy density and the energy 
are finite. At the origin the gauge and Higgs field functions satisfy the boundary conditions 

f A (0) - 1 = f B (0) = H'(0) = K(0) = . (41) 

At infinity the gauge and Higgs field functions lie on a circle 

f A (oo) + if B (oo) = exp(^(oo)) , H{oo) + iK(oo) = exp(t^-) , (42) 

where 9(oo) is an unknown function of ^. Therefore we choose the boundary conditions 

/a(oo) = f' B (oo) = H'(oo) = K'(oc) = . (43) 
For the gauge function $ we choose the boundary conditions 

$(0) = , $'(oo) = , (44) 
consistent with the boundary conditions of the gauge and Higgs field functions. 

3.3 Results 

In the calculations presented we have employed the gradient formalism with constraint 
on the distance, since it appeared numerically far more stable than the explicit gradient 
approach. Let us now discuss the sphaleron barrier as obtained in the gradient approach. 
Barrier 
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We first compare the gradient path with the minimum energy path for small Higgs 
boson masses. In Fig. 1 we show the energy as a function of the Chern-Simons number for 
Mh = Mw for the gradient path and the minimum energy path. The minimum energy 
path barrier is steeper with respect to the Chern-Simons number than the gradient path 
barrier. This picture reverts, when we consider the energy as a function of the pathlength 
/, defined by 

rf r(bi) sphaleron 

1 = df- df . (45) 

J vacuum J vacuum 

Note, that the pathlength / is shifted such that the sphaleron or bisphaleron has 1 = 0. In 
Fig. 2 we show the energy as a function of the pathparameter I for Mh = My/ along the 
gradient path and the minimum energy path. With respect to the pathparameter / the 
gradient barrier is steeper. 

In Fig. 3 we present two configurations with the same Chern-Simons number, one along 
the gradient path and one along the minimum energy path. For the latter configuration 
the asymptotic values of the functions are already closer to the vacuum values. 

In Fig. 4 we show the energy, obtained with the gradient method, as a function of 
the Chern-Simons number for Mh = 15My/. Now there are three extrema of the energy 
functional, the two degenerate bisphalerons and the symmetric sphaleron. The right 
bisphaleron barrier is obtained from the left one by the transformation Nqs — > 1 — iVcs 
and E — > E. All three barriers, the lower asymmetric bisphaleron barriers and the higher 
symmetric sphaleron barrier, are smooth in the gradient approach, in contrast with the 
bifurcations encountered along the extremal energy path. Note, that the asymmetric 
bisphaleron barrier culminates at the bisphaleron and has on one side a very steep fall-off 
|f?5|] . In Fig. 5 we show these energy barriers as functions of the pathparameter I. 

The effective mass m is shown in Fig. 6 as a function of the Chern-Simons number for 
Mh = My/ and Mh = 15My/. For small Higgs boson masses, we find a smooth effective 
mass, qualitatively similar to |50| . But for large Higgs boson masses the effective mass 
develops a sharp peak, when the Higgs field crosses zero at spatial origin. This point 
coincides with the symmetric sphaleron, but not with the bisphaleron. In the vicinity of 
this peak the potential energy falls off steeply. This steep fall-off occurs only on one side 
(the side of the peak) of the bisphaleron, but on both sides of the symmetric sphaleron. 

By computing the distance of the (symmetric) sphaleron to the vacuum we obtain 
an estimate of how good the paths are. For Mh = My/ the distance is 40, while the 
pathlength of the gradient path is 46, and the pathlength of the minimum energy path is 
50. Thus the minimum energy path has a longer pathlength than the gradient path. 

Tunneling amplitude 

As a related criterion for the quality of a path let us now consider the associated 
semiclassical tunneling amplitude, determined by exp(— R ) |28], |29], [30], [JT], Q 

R = f b d\J2m(X)(V(X) - E) , (46) 

J a 

where V(X) is the potential energy and E is the energy of the classical turning points a 
and b [BO, BT1. (Note that the integral is independent of a reparametrization.) Considering 
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vacuum to vacuum transitions, (with respect to the above metric, Eq. (p7|),) the exponent 
Ro of the tunneling amplitude is a line integral along the path in configuration space, 




^o= / df\l-^, (47) 
where 

df=\imd(fj) , (48) 
/-/ 

V(X) = Eb, and E — 0. Employing the radial gauge in the calculations, the tunneling 
amplitude is determined by Rq as given in Eq. (f|7|) with d(f,f), Eq. (p7|), replaced by 
d(f, f), Eq. (|35|) , where we now interpret the sequence of configurations along the path as 
gauge transforms of those configurations for which d(f, f) equals d(f, f). 

For M H = M w we find i? = 1-575 and 1.836, in units of for the gradient path and 
the minimum energy path, respectively. For Mh = l5Mw we find Ro = 1.525 (1.533), 
in units of for the bisphaleron barrier (symmetric sphaleron barrier) in the gradient 
approach. In Fig. 7 we show Ro as a function of the Higgs boson mass in the gradient 
approach. We observe, that R ~ 1.5, in units of 8Tr 2 /g 2 , fairly independent of the Higgs 
boson mass, and the bisphaleron transition rate is slightly higher than the sphaleron 
transition rate. 



4 Level Crossing 

Let us now consider the fermionic level crossing phenomenon along the sphaleron barrier. 
We study the fermion mode in the background field of the barrier as well as selfconsistently, 
and compare the gradient path to the extremal energy path. 

4.1 Energy Functional 

To retain spherical symmetry we consider only fermion doublets degenerate in mass. The 
corresponding spherically symmetric ansatz for the fermion eigenstates is the hedgehog 
ansatz, 

q L (f,t) = e~ iwt Ml[G L (r)+ia-rF L (r)]xh , (49) 
q R (f,t) = e-^Ml[G R (r) - id ■ fF R (r)] Xh , (50) 
where the normalized hedgehog spinor Xh satisfies the spin-isospin relation 

&Xh + TXh = . (51) 
Under the residual gauge transformation, Eq. [15], the fermion functions transform as 

F L + iG L — ► exp(i®)(F L + zG L ) , 

F R + iG R — > F R + iG R . (52) 
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The fermionic energy functional reads 
Et = 
- 2 



AttM w / dxx 
Jo 

1-fA 



2 2 
FrGr — G R F R H — F R G R + F L G L — G L F L H — F L G L 

a; x 



x x " 2x 



G L F L + i±{Gi - Ft) + ^(Gi + i^) 



+ 2M F H(G R G L - F R F L ) - 2M F K{F R G L + FlGr' 
where the fermion mass is expressed in units of M\y 

Mr = M F /M W . 



(53) 



(54) 



The fermion functions need to be normalized. When N fermions occupy the eigenstate 
the normalization condition is 



POO 

4tt / dxx 2 (G 2 R + F R + G 2 L + F?) = N 
Jo 



(55) 



4.2 Background Field Calculation 

Let us first consider the fermions in the background field of the sphaleron barrier. We find 
the set of coupled equations 



23| 



CjG l - F' L 



-F L + - — —Fl - —Gl + M F (-HG R + KF R ) = 

iy nn >~Y* 

iAj ju ju 



uF L + G' L + ^—^G L + f -^F L + M F (HF R + KG R ) = 







(56) 

(57) 

(58) 
(59) 

(60) 
(61) 

G R {0) = c R , G L (0) = c L , (62) 

where c R and c F are unknown constants, subject to the normalization condition fl5"5"|). At 
infinity all fermion functions vanish 



- uG R + F' R + -F R + M F (HG L - KF L ) 

CuF R + G' R + M F {HF L + KG L ) = 
where u is the fermion eigenvalue u in units of Mw 



Mi 



w 



At the origin the fermion functions satisfy the boundary conditions 

Fr(0) = F L (0) = , 



and 



F R (oo) = F L (oo) = G R (oo) = G L (oo) = . 



(63) 
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Let us now consider the case M# = My/. The fermion eigenvalue along the gradient 
path is shown in Fig. 8 as a function of the Chern-Simons number for the fermion masses 
Mp = lOMw, Mp = Mw and Mp = Miy/10. For comparison the fermion eigenvalue 
along the minimum energy path is also shown |T3 . For small fermion masses the fermions 



are bound only in the vicinity of the sphaleron. Here we find qualitatively the same 
behaviour of the fermion eigenvalue. For heavier fermions the eigenmode reaches the 
continua later along the gradient path than along the minimum energy path. 

For large values of the Higgs boson mass, when the barrier culminates at the bi- 
sphaleron, the fermion eigenvalue in the gradient approach is a monotonic function of the 
Chern-Simons number as shown in Fig. 9 for Mh = 15My/ and Mp = lQMyy, Mp = 
My/ and Mp = M^/IO. This is in contrast with the extremal energy path, where the 
bifurcations along the path also lead to bifurcations of the fermion eigenvalue [pO~ . Note, 



that the eigenvalue along the right bisphaleron barrier is obtained by the transformation 
^cs — > 1 — Nqs and uj — > —uj. 

In Fig. 10 we show the dependence of the zero eigenvalue of the fermions on the Chern- 
Simons number and on the fermion mass for the gradient path and the extremal energy 
path for Mn — 15Mjy. Depending on the Higgs boson mass, the zero mode approaches 
a limiting value for large fermion masses in the gradient approach. In contrast, along the 
extremal energy path the zero mode occurs for large fermion masses only at the sphaleron, 
i. e. at iVcs = 1/2. For small fermion masses the level crossing occurs for both methods 



in the vicinity of N cs = 1/2 Ml, M. 



4.3 Selfconsistent Calculation 

Let us now study the gradient path over the sphaleron barrier in the presence of fermions. 
We proceed analogously to our previous calculation p3f , but compute the barrier with the 
gradient method. We arrive at the same set of equations for the fermion fields, while we 
have to add the source terms 

+ g 2 xF L G L , (64) 
+ \M&l - Fl) , (65) 

+ g -^ L {G R G L -F R F L ) , (66) 
q 2 M f 

- g -^L(F R G L + F L G R ) , (67) 

to the right hand side of the boson field equations for f^, H, and K, respectively. 

Let us first consider small Higgs boson masses, where only the sphaleron barrier exists. 
As before, when studying the fermion eigenmode along the minimum energy path |23 



we observe that fermions with small masses have little influence on the shape of the 
barrier, while heavy fermions deform the barrier considerably. However, for very large 
fermion masses, the bifurcations, which we observed previously along the extremal energy 
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path, are no longer present along the gradient path. In the gradient approach the barrier 
decreases monotonically to both sides of the sphaleron, as shown in Fig. 11 for Mh = Mw 
and M F = 75 M w . 

As before [E3] we observe, that the fermion eigenvalue deviates litte from the eigenvalue 
of the background field calculation for small fermion masses, also for heavier fermions the 
path does not lead to a free fermion solution but to a bound state, a nontopological soliton. 
The selfconsistent fermion eigenvalue along the gradient barrier is shown in Fig. 12 for a 
heavy fermion with Mp = lOMw for the Higgs boson mass M# = Mw, and compared 
to the eigenvalue of the background field calculation. In the selfconsistent calculation the 
soliton is approached for Ncs — > 0. 

Let us now turn to large values of the Higgs boson mass, where we expect two bi- 
sphaleron barriers beside the sphaleron barrier. The presence of the fermions lifts the 
degeneracy of the two bisphalerons for finite fermion masses. Considering now the total 
energy, consisting of the bosonic energy and the fermion eigenvalue as encountered along 



the path over the barrier, we expect |3q| , that the energy of the left bisphaleron first 
increases as E = + Mp, while the energy of the right bisphaleron decreases as E = 
Ei, — Mp. (The left sphaleron is encountered along the barrier before the level crossing, 
i. e. the fermion is still in the positive continuum, while the right sphaleron is encountered 
after the level crossing, i. e. the fermion is in the negative continuum [pOfl .) At a critical 
value of the fermion mass the fermion becomes bound, Mf w A^J Mh / Mw — 12 GeV [j37jj . 
Then the energy of the left bisphaleron decreases, while the energy of the right bisphaleron 
increases. Interestingly, a bifurcation occurs at a moderate value of the fermion mass. 
The right bisphaleron merges with the sphaleron at a critical value, beyond which only 
the left bisphaleron solution exists. This curious feature of the selfconsistent treatment is 
demonstrated in Fig. 13 for Mh = 15Mw (for one bound fermion). 

The selfconsistent fermion eigenvalue for a heavy fermion with Mp = lOMw for the 
Higgs boson mass Mh = 15Mw along the left gradient bisphaleron barrier is shown in 
Fig. 12. Note, that this is the only selfconsistent barrier for this fermion mass. Here, the 
vacuum is approached for iVcs — * 0, because the soliton exists only for higher fermion 
masses p3fl. 



5 Conclusion 

We have applied the gradient approach to obtain the sphaleron barrier. The gradient 
approach produces a path of steepest descent with respect to a given metric. A natural 
metric on the space of field configurations is implied by the kinetic energy term of the 
Lagrangian []27| . We have formulated this metric in a gauge invariant way [[34[ . 



We have presented the formalism of the explicit gradient approach and of the gradient 
approach with a constraint on the distance. For technical reasons we have used the latter 
approach in the numerical calculations. 

Since the bifurcations along the extremal energy path have largely motivated this 
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study, we have compared the sphaleron barrier obtained along the gradient path to the 
one along the extremal energy path. For small values of the Higgs boson mass, there is 
only one sphaleron solution, the symmetric sphaleron 0, |l|, [14|. Here both approaches 
lead to a smooth barrier. The gradient path barrier is steeper with respect to the path- 
length /, defined via the metric, while the minimum energy path is steeper with respect 
to the Chern-Simons number. But the semiclassically calculated tunneling amplitude, 
~ exp(— Rq), is bigger for the gradient path, e. g. for M# = M\y R is smaller by 10% for 
the gradient path. 

For Higgs boson masses larger than ~ 1 TeV new asymmetric sphaleron solutions with 
lower energy appear, the bisphalerons, |13|, 14] . When bisphalerons exist, the extremal 



energy path has bifurcations and culminates not at the bisphaleron but at the symmetric 
sphaleron. In contrast, the gradient approach leads to smooth barriers, a lower asymmetric 
bisphaleron barrier and a higher symmetric sphaleron barrier. The asymmetric barrier has 
a steep fall-off on one side. This fall-off is related to a peak in the effective mass, when 
the Higgs field passes zero (at the origin). The semiclassical tunneling amplitude is fairly 
independent of the Higgs boson mass, but slightly bigger along the bisphaleron barrier 
than along the sphaleron barrier. 

To exhibit the level crossing phenomenon we have calculated the valence fermion mode 
along the gradient approach barriers. Since the barriers are smooth in the gradient ap- 
proach, also the fermion eigenvalue along the barriers is smooth. This is in contrast to the 
extremal energy path, where the bifurcations of the barriers were reflected in bifurcations 
of the fermion eigenvalue |[20|| . The fermion eigenvalue decreases monotonically from the 
positive continuum to the negative continuum along the gradient path in the background 
field approximation, even for large values of the Higgs boson mass. 

When fermions are coupled selfconsistently to the boson fields the fermion mass is of 
importance. For small fermions masses there is hardly any change with respect to the 
background field calculations, while for heavy fermions the barriers deform considerably. 
Notably, for large Higgs boson masses, two of the three barriers, the sphaleron barrier 
and the right bisphaleron barrier, merge and disappear already for moderate values of the 
fermion mass, leaving as the only barrier the left bisphaleron barrier. 
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6 Figure Captions 



Figure 1: The total energy (in TeV) is shown as a function of the Chern-Simons number 
iV cs along the gradient path (solid) and along the extremal path (dotted) for M H = M w . 

Figure 2: The total energy (in TeV) is shown as a function of the pathparameter / 
along the gradient path (solid) and along the extremal path (dotted) for Mh = Mw- 

Figure 3: The gauge field functions J'a and fs are shown with respect to the dimen- 
sionless variable x for the configurations with Nqs = 1/2 along the gradient path (solid) 
and along the extremal path (dotted) for M H = M w . 

Figure 4: The total energy (in TeV) is shown as a function of the Chern-Simons num- 
ber iVcs along the symmetric sphaleron path (solid) and along the asymmetric bisphaleron 
path (dotted) in the gradient approach for Mh = IbMw- 

Figure 5: The total energy (in TeV) is shown as a function of the pathparameter I 
along the symmetric sphaleron path (solid) and along the asymmetric bisphaleron path 
(dotted) in the gradient approach for M H = 15M W . 

Figure 6: The effective mass (in units of 1/Mw) is shown as a function of the Chern- 
Simons number Ncs i n the gradient approach along the symmetric sphaleron path for 
Mh = Mw (solid) and Mh = l5Mw (dot-dashed), and along the asymmetric bisphaleron 
path for M H = 15M W (dotted). 

Figure 7: The transition amplitude -R (i n units of 8ir 2 /g 2 ) is shown as a function of 
the Higgs boson mass Mh (in units of Mw) in the gradient approach for the symmetric 
sphaleron path (solid) and for the asymmetric bisphaleron path (dotted). 

Figure 8: The fermion eigenvalue (in units of M F ) is shown as a function of the 
Chern-Simons number Ncs in the background of the sphaleron barrier for the fermion 
masses M F = M^/IO, M F = M w and M F = 10M W along the gradient path (solid) and 
along the extremal path (dotted) for Mh = Mw- 

Figure 9: The fermion eigenvalue (in units of M F ) is shown as a function of the 
Chern-Simons number Ncs in the background of the bisphaleron barrier for the fermion 
masses M F = Mvk/10, M F = M w and M F = 10M W along the gradient path (solid) for 
M H = 15M W . 

Figure 10: For the fermion zero mode the dependence of the fermion mass (in units 
of Mw) on the Chern-Simons number Ncs is shown along the gradient path (solid) and 
along the extremal path (dotted) for Mh = 15Mw- 

Figure 11: The energy (in TeV) is shown as a function of the Chern-Simons number 
N C s f° r the fermion mass M F = 75M W in the selfconsistent calculation along the gradient 
path (solid) and along the extremal path (dotted) for M H = M w . 

Figure 12: The fermion eigenvalue (in units of M F ) is shown as a function of the 
Chern-Simons number Ncs in the gradient approach along the sphaleron barrier for Mh = 
M w and along the bisphaleron barrier for M H = 15M W for the fermion mass M F = 10M W 
in the background field calculation (dotted) and in the selfconsistent calculation (solid). 

Figure 13: The energy (in TeV), including the fermion eigenvalue, is shown as a 
function of the fermion mass M F (in units of M w ) in the selfconsistent calculation for 



19 



the left bisphaleron (dotted), for the right bisphaleron (solid), and for the sphaleron 
dashed) for M H = 15M W . 
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